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Abstract 



Classical relativistic field theory is applied to perfect and magneto- 
hydrodynamic flows. The fields for Hamilton's principle are shown 
to be the Lagrangian coordinates of the fluid elements, which are 
potentials for the matter current 4-vector and the electromagnetic 
field 2-form. The energy momentum tensor and equations of motion 
are derived from the fields. In this way the theory of continua is shown 
to have the same form as other field theories, such as electromagnetism 
and general relativity. 

Waves are treated as an example of the power of field theoretic 
methods. The average or background flow and the waves are con- 
sidered as two interacting components of the system. The wave- 
background interaction involves the transfer of energy and momen- 
tum between the waves and the average flow, but the total energy 
and momentum are conserved. The average Lagrangian for the total 
wave-background system is found by expanding the Lagrangian about 
the background flow and averaging over the phase. The total energy- 
momentum tensor is constructed, and the conservation of energy and 
momentum are discussed. Varying the wave amplitude gives the dis- 
persion and polarization relations for the waves, and varying the phase 
gives the rays and conservation of wave quanta (or wave action). The 
wave quanta move with the group velocity along the bi-characteristic 
rays. 

Part I 

Field Theory for Fluids 

1 Introduction 

The theory is both Lorentz and generally covariant on space time M, so the 
physical and geometric quantities are independent of the coordinate system, 
i.e. they are tensor fields. We use coordinate free notation, and units so that 
the speed of light = 1. The natural mathematical framework is differential 
geometry 00E!- 

The theory is introduced in Part 1 and is applied to waves in Part 2. 
Relativistic thermodynamics is discussed in §2. In §3 the Lagrangian coordi- 
nates are shown to be the three scalar fields for Hamilton's principle. They 
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are also the potentials for the matter current 4-vector and the electromag- 
netic field 2-form in magnetohydrodynamics (MHD) and plasmas. In §4, 5, 
and 6, the energy momentum tensor and the conservation equations for rest 
mass (or the number of particles) and for the energy and momentum, are 
discussed for perfect fluids and MHD. In Part 2, the theory is applied to 
waves. The fields for the waves are their Lagrangian coordinates, which are 
the oscillation centers. 

2 Relativistic Thermodynamics of an Ideal 
Gas 

A proper quantity is one that moves with the flow (i.e. is measured with 
respect to a comoving coordinate system) and so is coordinate independent. 
It is a scalar field, which is a rank tensor field (tensors include scalars and 
vectors). I use coordinate independent notation to emphasize the fact that 
tensors are independent of the coordinate system. The following thermody- 
namic quantities are scalar fields since they are all proper quantities. 

• n' — proper number density 

• m = rest mass of fluid particles 

• n = n'm = proper density of rest mass 

• ne — internal (thermal) energy density 

• e = thermal energy per unit mass 

• s = specific entropy 

• p = n + ne = total energy density 

• T = temperature 

• p = nT/m = isotropic pressure 

• 7 = c p /c v = d\np/d\nn\ s = the ratio of specific heats 

• / = (p + p)/ n — l + e + p/n = the index or relativistic specific enthalpy 
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If n and s are the independent thermodynamic variables, the first law of 
thermodynamics is the conservation of energy: 



Tds = de + pd(^j (1) 

If the flow is isentropic, ds = and s is constant. Then the thermodynamic 
quantities are functions of n alone, and equation (1) gives 

P = n 2 ^- dp = ndf ^ = f. (2) 

dn dn 



3 The Mathematical Formalism 

3.1 Three Coordinate Systems 

Assume that there is no gravity so that spacetime M is flat. Consider the 
following three coordinate systems for M: 

1. The observer's global orthonormal coordinate system 

2. A local comoving coordinate system (CMCS) y = (y° , y l , y 2 , y 3 ) at 
each point x G M. This coordinate system is also orthonormal, with 
u = d/dy°. 

3. A global CMCS z = (z°, z 1 , z 2 , z 3 ). z is not orthonormal, so the metric 
g is not diagonal. In this coordinate system, u \\ d/dz°, and z° is 
closely related to the proper time for a fluid element. The z k are the 
fields which label the worldlines. 



3.2 Differential p-Forms 

Let A P (M)(A P (M)) be the set of p-forms (p- vectors) on M. A contravariant 
vector field v G Ai(M) or a 1-form a G A X (M) can be written in terms of 
the dual coordinate basis e a = d/dx a and e a = dx a as 

v = v^e^ a = 

The exterior derivative d: A p — > A p+1 of a is da = V A a. A generalized 
inner product on forms and p-vectors can be defined as the adjoint of the 



5 



exterior product: 



( | ) p : A q x A p -> A q - p : (a, (3) -> (a|/3) p = a(/3) g > p, denned by 
(a(/3)| 7 ) = a(/3,7) = (a|/3 A 7) where a G A 9 , (3 G A p , 7 G A q ~ p . 

The Hodge operator * is defined by * : A p — > A 4_p : /3 — > */3 

a A */3 = (a|/3)r ^ */? = (r\(3) p = t((3), where a, (3 G A p , (3) 
r = *1 = a/- det g e, and e = dx° A dx 1 A da; 2 A da; 3 . 

r is the volume element 4-form, e is the 4D Levi-Civita tensor, and g is the 
metric. If a; is a orthonormal coordinate system, detg = — 1, and r = e = d 4 a;. 

3.3 The Scalar Fields 
3.3.1 Equation of continuity 

Compare the problem of a fluid in 2D spacetime (t, x) = (x°, x 1 ) with electro- 
magnetism (EM). (The following is also true in 4D spacetime). The equation 
of continuity dj = is analogous to Maxwell's homogeneous equation dF = 0. 
j is the dual of the matter current vector J = nu = *j, and j = * J = * * j. 
We have 

V • J = *{dj) = dj = j = dz (4) 

Thus the equation of continuity implies the existence of the scalar field z(t, x) 
such that j = dz, just as in EM where we have the vector field A and F = dA. 
In both cases we use the fact that a closed form is exact. In 4D spacetime, 
there are 3 scalar fields, and 

dj = j = dz 1 A dz 2 A d^ 3 

Often a density r is introduced so that z % can be a Cartesian coordinate 
system, and then j = r dz 1 A dz 2 A dz 3 . The invariant density is 

n = y/UU) 

Consider a fluid in flat spacetime M C IR 4 (no gravity). We need to 
express the Lagrangian in terms of the fields and their first partial derivatives. 
For the initial value problem, x° = t > 0. Define the initial configuration on 
the observer's 3D hypersurface SCR 3 defined by x° = t = 0. The equations 
of motion determine the flow on M for t > 0. 
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The three scalar fields are the Lagrangian coordinates Z = (z 1 ,z 2 ,z i ). 
They are a coordinate system on £ and move with the fluid, and are funda- 
mental to the theory of all continua. The flow is determined by the projection 
Z of a fluid element at x E M onto its initial position Z(x) on £: 



and the unit 4-velocity u is tangent to the worldlines of the fluid elements 
(i.e., the streamlines). The z k have the following properties: 

1. They are potentials for the matter current 4- vector J = nu and the 
2-form of the electromagnetic field F. 

2. They are scalar fields which label fluid elements and so are constant 
along worldlines (dz k (u) = 0). 

3. They label worldlines, so varying the z k varies the worldlines of the 
fluid elements. 

4. By using the fields z k we automatically satisfy the equation of conti- 
nuity dj = and conserve matter. 

3.4 The Matter Current and Proper Number Density 

Let n : £ — > ffi: z k — > n (z k ) be the number density on the initial hypersur- 
face E. The 3- form of a fluid element on the £ is 



is the total number of particles (or the total rest mass). The flow is described 
by the matter current 4- vector J = nu E A^M), or its dual 3-form j = *J E 
A 3 (M). j = Z*(j ) is the pull-back of j from E to M by Z. 



Z: M ^ E: x -> Z(x) = (z\x), z 2 (x), z 3 (x)) where 
z k (0, x l ) = x k and dz k (u) = u a — = 




jo = n dz 1 A dz 2 A dz 3 and = j 



E 



Z*: A 3 (E)^A 3 (M): jo^^*(jo) 




.a. 



',/3-y 



where e a/37 = e a A A e 7 and r = no o Z. 
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The matter current 4-vector is 

J = *j = e(j) = e a ^ s j a/3l e delta /3\ = J s e s . (5) 
Thus the components of J and j are 

. d(z\z 2 ,z 3 ) . . d{z\z 2 ,z 3 ) 1 

J = J123 = r = 7^ J = -J023 = ^ a/ n 2 3 x = 7™> 

2 . d^ 1 ,* 2 ,* 3 ) 2 3 . <9(^W 3 ) 3 

J = 3013 = r ^wj = 7m; J = _J012 = " r ^W) = 7 ™ 

These expressions show that the fields z % are indeed the potentials for the 
components of J. They also give the 3- velocity v l = J l /J° in terms of the 
z\ and the components of the 4-velocity u; u° = 7 = l/sqrt(l — v 2 ) and 
u l = 7f l . 

The proper number density n is 

n 2 = = (*J| * J) = r 2 ^ 1 A dz 2 A dz^dz 1 A d^ 2 A dz 3 ) (6) 

= r 2 det(d^|d^') 

r is constant along streamlines since the equation of continuity is 

dj = d(*J) = *(V • J) = dr A d 3 z = dr(u) = r = r(z i ) (7) 

The equation of continuity follows from djo = 0, since there are no 4- forms 
in S: 

dj = d{Z*j Q ) = Z*(dj ) = 0, since d o Z* = Z*o d 

A current J with V • J = is said to be conserved. Associated with any con- 
served current is a conserved 'charge', which is the total number of particles 
(or total rest mass) 

N — ( j — [ jn I J°d 3 x 



J — Jo 
s Jt, jt, 



3.5 Decomposing Vectors and Tensors 

The metric g and 4-velocity u define a spatial projection operator P = g + uu 
such that P(u) = 0, P = P T , P 2 = P, Tr(P) = 3, and P = diag(0, 1, 1, 1) in 
a local CMCS. P projects a vector field onto the 3D space orthogonal to u at 
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any x G M. A vector field A can be written as the sum of a vector parallel 
to u and a vector orthogonal to u: 

A = -A{u)u + P(A) = au + A'. 

a is the temporal part of A and A' is the spatial part, a and A' are proper, 
coordinate independent quantities, a is a scalar field and A' is a spatial vector 
field. 

If T is a second order tensor field, decompose it as follows; 

T = T(u, u) - uP{T{u)) - P{T\u))u + P-T-P 
If T is the energy-momentum tensor, 

• T(u, u) = proper energy density 

• —P(T(u)) = proper energy flux 

• —P{T T {u)) = proper momentum density, and 

• PTP = proper stress tensor 

Only proper quantities appear in the equations of motion, but the observer 
measures quantities with respect to his coordinate system. Observed and 
proper quantities are related by a Lorentz transformation with the 3- velocity 
v(x) of the fluid in the observer's x coordinate system. 

4 Energy-Momentum Tensor and 
Conservation Equations 

We must write the Lagrangian L in terms of the fields z k and their first 
partial derivatives dz k , and calculate the energy-momentum tensor (EMT) 

T = Lg-DL where D = dz k ® ^ = e v ; = *jA ( 8 ) 

The fields are scalars so T = T T and DL are symmetric. DL is spatial since 
dz k (u) = so 

DL(u) = (DLy(u) = dz k (u)— = 
9 



The proper energy density is T(u, u) = —L. The proper stress tensor is 
PTP = LP — DL. The proper energy flux P(T(u)) = = proper momentum 
density. This will not be true when waves are introduced, as things become 
more complex. 

The Euler-Lagrange equations which follow from varying the fields z k (i.e. 
varying the worldlines of the fluid elements) are 

9L „ T , dL 

- V • L k = where L k = ^—re^ 
oz" ozf, 

The Lagrangian has no explicit dependence on x, so the Euler-Lagrange 
equation is equivalent to conservation of energy and momentum: 

V • T = v^t;^ = (9) 
Decomposing equation Q with respect to u gives 

u ■ (V • T) = and P(V • T) = (10) 
These are the equations for conservation of energy and momentum. 



5 Perfect Fluids 

The Lagrangian for a perfect fluid is 

L = -p = -n{\ + e) (11) 

Equations © and (jHJ) for n and T yield Dn = nP and the well known EMT 
for a perfect fluid: 

T = Lg - DL = puu+ pP (12) 

The proper energy density is T(u,u) = p. The stress tensor is isotropic; 
PTP = pP, so T is diagonal in a local CMCS; T = diag(p,p, p,p). The 
equations of motion are 

V ■ T = nfa + dp + itV • (nfu) = 

where a = Viz (it) is the acceleration. 
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5.1 Conservation of Energy 

u ■ (V ■ T) = dp(u) - ndf(u) = -nTds(u) = (13) 

since V ■ (nu) = 0. This is the thermodynamic equation (1) for adiabatic 
flow (ds(u) = 0). Any two of the following three imply the third: 

V ■ (nu) = u ■ (V ■ T) = ds(u) = 

For isentropic flow (ds = 0), equation (|T3*j) for the energy is dp = ndf . 

5.2 Conservation of Momentum 

P(V ■ T) = nfa + P(dp) = (14) 
For isentropic flow, dp = ndf, and the equation of motion (|T4*|) is 

a = -P(df)/f (15) 

6 Magnet ohydro dynamics (MHD) 
6.1 The Electromagnetic Field F 

Decomposing the electromagnetic (EM) field 2-form F with respect to u yields 

F = u A E + *(u A B) 

where E = —F(u) is the proper electric field and B = *F(u) is the proper 
magnetic field. Both E and B are spatial vector fields, since F is anti- 
symmetric; E(u) = -F(u,u) = and B(u) = *F(u,u) = 0. A MHD 
fluid has infinite conductivity which implies a frozen in magnetic field, i.e. 
E = -F(u) = 0. Thus F = PFP is spatial, and in a global CMCS z, 

F = *(u A B) = e a b c b a dz b A dz c /2 (16) 

where e a bc is the 3D Levi-Civita tensor and b a are the components of F in 
the z coordinate system. So we have F in terms of the fields z k . Maxwell's 
homogeneous equation 

dF = db a (u) = 

So the 6 a (z 1 , z 2 , z 3 ) are constant on streamlines (the field is "frozen in"). 
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6.2 The Lagrangian and the Fields 

The Lagrangian for the magnetic field in MHD is L b = — (F\F)/2 = —B 2 /2. 
Equation (fH)|) gives F and thus L b in terms of the fields z 1 and their first 
derivatives dz l \ 

U = -e abc e def b a b d (dz b \dz e )(dz c \dz f )/4 (17) 

The total Lagrangian for MHD is the sum of L b and the Lagrangian for the 
matter (fTTj) L m = —p: 

L = L m + L b = -p- B 2 /2 (18) 

6.3 The Energy Momentum Tensor (EMT) 

From (HHJ) and © we find that DL b = BB - B 2 P, so the EMT for the 
magnetic field is 

T b = L b g - DL b = B 2 (uu + P) - BB 
Using (H2J) for T m gives the total EMT: 

T = T m + T b = (p + B 2 /2)uu + (p + B 2 /2)P - BB (19) 
If ei || B in a CMCS, T is diagonal: 

T 00 = p + 5 2 /2 T u =p- B 2 /2 T 22 = T 33 = p + B 2 /2 

6.4 The MHD equations 

The equations for MHD are V • (nu) = and 

V • T = V • [(n/ + B 2 )u}u + (n/ + S 2 )a + Vp + (Vfi) f (fi) (20) 
-(V • 5)5 - V5(.B) = 
dF = O V • *F = V • (Bu - uB) = (V • B)« + Vn(5) 

-(V • u)S - VB(u) = (21) 

6.4.1 Conservation of Energy 

The energy equation for MHD is the projection of equation (J2*U|) onto u; 
u ■ (V • T) = 0. It is the same as the energy equation (J13)) for a perfect fluid. 
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6.4.2 Conservation of Momentum 

The MHD equation of motion is the spatial part of (f2*Uj) : 

P(V ■ T) = (nf + B 2 )a + nP(df) - (V ■ B)B - P(dB(B)) = (22) 

Maxwell's inhomogeneous equation relates the electric current 4- vector K to 
B = *F(u): 

K = *d*F = *d(B A u) = *dB(u) - *du(B). 
The MHD equation of motion can be written as 

P(V • T) = (nf + B 2 )a + nP(df) - K x B = 
where is Lorentz force is K x B = P(K) x B = *(u A K A B). 

Part II 

Waves in Continua 

7 Formulating the Problem 

The following theory of waves in continua gives the dispersion relation and 
polarization of the various modes, the bi-characteristic rays, the equation 
governing the wave amplitude, the ray equations which determine the wave 
vector along the rays, and the EMT of the waves, the average background 
flow, and interaction between the waves and the average background flow. 
There are three 4-vector fields associated with any wave mode. They are 

1. The wave 4-vector k = uju + k, where u = —k(u) and k = P(k) are the 
proper frequency and wavenumber vector. 

2. The phase 4- velocity = u+v p is the 4- velocity of a surface of constant 
phase. P{v<f) — v p — ujk/k 2 is the (proper) spatial phase velocity with 
respect to the fluid, i.e. as seen in a CMCS. 

3. The group 4- velocity g = u + v g is the 4- velocity of wavefronts (or of 
wave quanta or wave packets). P(g) — v g — duj/dn is the spatial group 
velocity in a CMCS. 
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Waves in continua are described by the characteristic hypersurfaces (with 
normal k) and bi-characteristic rays (integral curves of w). k satisfies the 
dispersion relation. Wavefronts can be thought of as discontinuities in the 
solutions to the equations of motion. Consider plane waves with surfaces of 
constant phase given by (p(x) = (k\x) = k^x 1 *, so k = dip is the normal to 
the hypersurfaces (p(x) = constant. 

7.1 The Two Sets of Independent Coordinates and Fields 

Make the usual assumption that the average background flow varies little 
over a period or a wavelength of the waves (the WKB approximation). The 
following theory of the waves and their interaction with the background flow 
requires two sets of independent coordinates and fields, i.e. two maps from 
Euclidean to Lagrangian coordinates, together with a Lorentz transformation 
from the local CMCS y to the x coordinate system (CS) of the observer. The 
Lagrangian is a scalar field so it can be calculated using the coordinates that 
are most convenient. The two sets of coordinates and fields must be clearly 
defined. They are: 

1. z % : x — > z l (x) ; (x, z) refers to the average flow ; 

2. w l : y — > w l {y) — y l + (f> l (y) where (f> l {y) = a 1 cosip(y) ; 

(x, z) defines the slowly varying flow of the average position of the fluid 
elements (i.e. the flow of the oscillation centers), called the average or back- 
ground flow, x is the position of an oscillation center, rather than a fluid 
element. z l {x) is the initial position of the oscillation center of the fluid 
element whose oscillation center is at x. 

(y, w) defines the waves in a local CMCS y in the neighborhood of any 
x G M. The y CS is moving with the 4-velocity u of the fluid's oscillation 
centers, so that e = d/dy° — u in this CS. The fields w l (y) are the oscillation 
centers, the average position of the fluid element at y. <fi l (y) = a l cos(p(y) 
is the displacement due to waves. The spatial vector a = (a 1 , a 2 , a 3 ) is the 
wave amplitude. For a plane wave, the phase tp(y) = {k\y) is rapidly varying, 
but a and k = dip are slowly varying functions of x (not y). To treat the 
waves, we replace the fields w l by a % and ip, the appropriate fields for waves. 
Varying a gives the dispersion relation and polarization a = a/\a\. Varying 
ip gives the group velocity, the rays, and the equation for conservation of wave 
action (or wave quanta-dividing by Planck's constant gives a dimensionless 
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number that is conserved, which I call quanta, even though this is not a 
quantum theory). This determines how the wave amplitude evolves as the 
wave propagates through the fluid. 

8 The Average Lagrangian 

In the following, I will denote rapidly varying quantities with a prime, and 
averaged ones as unprimed. 

The (slowly varying) average Lagrangian L(z k ,dz k ,a,k;x) of the total 
wave-background system is derived as follows. Make a Taylor expansion 
about the background flow (of the oscillation centers) of the rapidly varying 
Lagrangian V of the actual flow, to order a 2 . Then average L' over the phase 
(p (Whitham 1974) to get the total average Lagrangian L = L + L w , where 
Lq is the Lagrangian of the background and L w is the wave Lagrangian. 
One can argue that the rapidly varying part of L would average out in the 
integration in Hamilton's principle anyway. 

8.1 Sound Waves in a Perfect Fluid 

Apply the coordinates (y, w) to (JEJ). We have 

n'(y) = r(w)\/det(dw i \dw : >) = n + 5n (23) 
dw % = dy l + d<p l — e l — a l k sin <p 

To order a 2 we have 

6n/n = —k ■ asimp — u 2 a 2 sm 2 ip/2 (24) 
Expand V = —pf = —n'[l + e(n')]about n' to order a 2 : 

r)T f) 2 T 

U = L (n + 5n) = L (n) + %^5n + %-^{5n) 2 /2 + • • • 

on an A 

where L (n) = —pin) = —nil + e) is the Lagrangian of the background, 

dL d 2 L fc 2 2 dp jT 

— — = -/, and — — = ; where c = — = — 

on on 1 n dp mj 

c is the speed of sound. 
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The average of any function f(ip) is 
1 f 2n 



27r Jo 



=► (0*> = O, <^>=y l 



(5n) = — nuj 2 a 2 /A, 



((Sn) 2 ) =n 2 (K-a) 2 /2 



The total average Lagrangian is 

L = (V) = L + L w = -nil + e) + nf[uj 2 a 2 - c 2 (k ■ a) 2 }/ A (25) 

For ideal gases, relativistic kinetic theory gives exact expressions for n, e, p, p, 7, 
and c as functions of T/m (Lightmann et.al. 1975). For a monatomic ideal 
gas, as T — > 0, 7 — » 5/3 and c 2 — > 5T/3m, and as T — > cxo,p — > 3p, so 
c 2 = rfp/rfp _> 1/3 and 7 -> 4/3. 

8.2 Magnetohydrodynamic (MHD) Waves 

Repeating this procedure for MHD, expand V = L' m + L' b = —p' — (F'\F')/2 
to order a 2 and average over <p to get 



The subscript refers to the background flow, w refers to the waves, m to 
the matter, and b to the magnetic field. 

Equation (j2*o} gives (L' m ) = L m = L 0m + L wm . The calculation for L b 
is analogous to that for (j23J). Using equation (fTo^) for F' = *(u A B') and 
substituting (y, w) for (x, z) in L' b gives 



is the part of L w due to the magnetic field. After averaging, set {w l ) = y l — > 
x l since x % is the position of the oscillation center of the fluid element. Thus 
the total average Lagrangian is 




(26) 



L' b = -e abc e def B a B d (dw b \dw e )(dw c \dw f )/A 
Putting dw l = e l — a l k sin ip in (|2*7j) and averaging yields 



(27) 



L fe = (L' b ) = -B 2 /2 + L bw , where 
4L bw = k 2 [(a ■ B) 2 - a 2 B 2 } + (B A a A k|S A a A /c) 



(28) 




(29) 
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8.2.1 Alfven Waves: B ■ a = k ■ a = 
For Alfven waves (|2*9j) gives 

2L W = (nf + B 2 )u 2 a 2 - (k ■ Bfa 2 (30) 

8.2.2 Magnetosonic Waves: B A a A k = 

For the fast and slow magnetosonic modes ()29|) gives 
4L W = [nf + B 2 -(a- B) 2 ]uj 2 a 2 - nfc 2 (K ■ a) 2 - n 2 a 2 [B 2 - (a ■ B) 2 } (31) 

9 The Dispersion Relation, Phase Velocity, 
and Polarization 

To discuss the dispersion and polarization relations write L w as a symmetric 
quadratic function of the wave amplitude a on the 3D space orthogonal to u 
(at each x G M): 

L w = C(a, a)/2 

where C is a 3 x 3 matrix. The fields for the waves are a and (p, the action- 
angle coordinates. The Euler-Lagrange equation for a is 

BL 

— ™ = c(a) = det C = 

oa 

which is the dispersion relation, det C = is a cubic for v 2 = u 2 /k 2 , so 
there are 3 solutions for v 2 and 3 corresponding solutions of C(a) = for 
the polarization a = a/\a\. The phase 4- velocity is = u ± v p . 

9.1 Sound Waves 

Writing equation (j23|) in the form L w = C(a,a)/2 yields 

C = nf{u 2 P - c 2 kk) 

P is the identity matrix of the 3D space orthogonal to u. If e\ || k, C is 
diagonal, the dispersion relation is detC = cj 4 (o; 2 — c 2 k 2 ) = 0. v p — 
when k ■ a = 0, corresponding to transverse modes which are comoving 
discontinuities and do not propagate. The longitudinal mode with k \\ a is a 
sound wave with 

to 2 = c 2 k 2 =^> v p = v g = ck (32) 
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9.2 MHD Waves 



Let ei || B, k in the (ei, e^) plane, and cos 9 = k ■ B. Then 





fC n 


C\2 





H 




C22 







\ o 





C33 



and det C = C 33 [CnC 2 2 — (Cu) 2 ]- The are given below. MHD waves 
fall into two classes, depending on the polarization a. When a 1 = a 2 = 0, 
a = is _L to the plane containing B and k. These are Alfven waves with 
dispersion relation C33 = O.When a 3 = and a lies in the plane of B and k, 
we have magnetosonic waves with dispersion relation 

C11C22 — (C12) 2 = 



9.2.1 Alfven Waves 

Alfven waves are transverse, with a _L both k and B, and equation (JHUj) for 
L w = C 33 a 2 /2 gives the dispersion relation 



v 



{B-kf _ , , fl . , 



I = f i P2 = w a cos2 ^ ! v p = cos $ /t where ^ A - 

nj + B z nj + B 

(33) 



9.2.2 Magnetosonic Waves 

5, «, and a all lie in the same plane, soBAnAa = (see (}3*Tj) for L w ). The 
dispersion relation for the fast (+) and slow(— ) (compressive) magnetosonic 
modes is det C — 0, which is a quadratic equation for the square of the phase 
velocity v± = uj 2 / k 2 : 

2C = {nf + B 2 )uj 2 P - u 2 BB - nfc 2 KK - k 2 (B 2 P - BB) 
2v 2 ± = c 2 + v 2 A - c 2 v 2 A sin 2 9 ± R where (34) 

R 2 = ( c 2 - V \ f + 2c 2 v 2 A sin 2 9\2 - c 2 - v\ + c 2 ^ sin 2 0/2] and 

v p = v ± k 

Note that t>_ < w + and c and f a appear symmetrically in v±. 
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9.3 Polarization 



The polarization of sound waves (a \\ k) and Alfven waves (a _L k) is clear. 
For the magnetosonic modes, let cos/5 = a ■ B. Then a± = cos/3ei + sin/3e2, 
and C(a) = =>- 

tan/3 =-^1 = -^ (35) 

The Lagrangian (|31|) and the dispersion relation (J34j) give the components 
of C as functions of 8, the angle between B and k. Then (J35|) gives the 
polarization (3(8) and «(#). The waves are discussed in detail in another 
paper (reference here). 



10 Wave Quanta, Group Velocity, and Rays 

Wave action can be replaced by the equivalent, more intuitive concept of 
wave quanta (let Planck's constant H — 1). Even though everything is con- 
tinuous, we will refer to the wave action as wave quanta. The discussion of 
wave quanta, group velocity, and the rays, is clarified by writing the wave 
Lagrangian as a symmetric function of k: 

L w = Q(k,k)/2 (36) 

where Q = Q Ms symmetric. The Euler-Lagrange equation for (p expresses 
conservation of wave quanta. Since <p has been averaged over, we have 

dip ax^ oo^ip ok 
where the wave quanta current 4-vector 
8L 

A = = -Q(k) = Nu + P(A) = N(u + v g ) = Ng (38) 

is a conserved current, 

N = -A(u) = ^ and P(A) = Nv g = (39) 

OUJ OK 

are the proper density and flux of wave quanta. The (proper) wave energy 
density E and group velocity are 
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Equation (J38|) for A shows that the wave quanta move with the group velocity. 

£ is the hypersurface defined by x° = t = 0, with normal dt = e° and 
volume element *(eo) = e 123 = dx 1 Adx 2 Adx 3 = d 3 x. The conserved quantity 
associated with the conserved current A is the total number of wave quanta 



10.1 Sound Waves 

Writing equation(j2SJ) as L w = Q(k, k)/2 shows that 



Q — nf(a 2 uu — c 2 aa)/2 

A = -Q(k) = N(u + v g ) = nfa 2 (uu + c 2 k) 

N = nfua 2 /2 = E/u 

v g = ck ; g = u± ck 



2Q = (nf + B 2 )a 2 uu - a 2 BB 
2A = -2Q(k) = (nf + B 2 )a 2 uju + a 2 (n ■ B)B 
N = -A(u) = (nf + B 2 )uja 2 /2 = E/u 
P(A) = Nv g = a 2 (K-B)B/2 

=>• v g — v aB and g = u± vaB 



N w : 




10.2 MHD Waves 



10.2.1 Alfven Waves 



Equation (|HU|) for L w gives 
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10.2.2 Magnetosonic Waves 

Define cos rj = k ■ &, and recall cos 9 = k- B and cos (5 = a ■ B, so 9 = (3 + 77, 
since the 3 vectors all lie in one plane. Equation (}3~T]) for L w gives 

(nf + B 2 sin 2 (5)a 2 uu — nfc 2 aa — B 2 a 2 sin 2 (3P 
-Q(k) = Nu + P(A) 
-A{u) = ua 2 {nf + B 2 sin 2 /3)/2 = E/u 
2iVt> g = n/c (k ■ a) a + -B a sin p/t 
(1 — t>i)c 2 cosi]±a± + v\ sin 2 /3±k 
v±(l — v\ cos 2 

v~ is confined to a narrow cone whose axis is B, but v£ can be in any 
direction, v^ 1 , like v±, is unchanged when c 2 and v\ are interchanged. Also 

2L W = Q(k, k) = k-A = Nk-g = 
vl = v g -v p ^v p < v g 

10.3 The Normal Cone, The Ray Cone, and the Ray 
Equations 

The bi-characteristic rays are the world lines of wave quanta and integral 
curves of the group 4- velocity g = u±v g or the wave quanta current A = Ng. 
There are two directions for the rays of each mode, corresponding to g — u+v g 
and g = u — v g . The 3D characteristic hypersurfaces {C : <f(x) = constant} 
are wavefronts, across which discontinuities occur. Since k = dip and k-g — 0, 
rays lie in characteristic surfaces. The wavefront for a disturbance at the 
origin x = divides spacetime M into two regions: behind or ahead of the 
wavefront. The wavefront can be thought of as the envelope of the surfaces 
orthogonal to the rays formed by a burst of wave quanta from the origin. 

Q defines two quadratic cones at each point of spacetime: the normal 
cone Q(k, k) = 0, which is generated by the rays which are integral curves 
of A, and the ray cone Q _1 (A,A) = 0, dual to the normal cone, which is 
generated by the wave-vector k. Each of the two cones is the envelope of 
the planes orthogonal to the rays of the other. It is useful to look at the 
spatial cross sections of these cones obtained by setting the time components 
equal to 1. In MHD, the surfaces have axial symmetry about the magnetic 
field. The cross section of the normal cone Q(k, k) — is called the normal 



2Q 
A 
N 

2P{A) 
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or wave- vector surface. The cross section of the ray cone Q~ 1 (A,A) = 
is called the ray surface. Its useful to plot v g {6) and v p {6). For a further 
discussion of the non-relativistic results see Courant and Hilbert Vol. II, §VI 
(1962). 

Since k = dip, dk = d 2 (p = 0. If k' = dp/dx 11 are the components of k in 
the x coordinate system (CS), 



dk 



dk' dk' 



dx^ 

at 



dx u 
duo 1 
dx 1 



dx 11 A dx u /2 = 



(41) 



and 



84 _ dK' 3 

dxi dx' 1 



where u' = uj'(k'(x),x) is the dispersion relation, k'^ = (— u', k^) are related 
to the proper components of k^ = (— lj, m) by a Lorentz transformation with 
the local fluid velocity v(x): 



uj' = j(uj + k ■ v) 



7(/C|| + ojv) ; and n' ± = k±_ (42) 



Consider the dispersion relation as a function of x, through both k and 
the inhomogeneities in the background flow: 



dL^(yk{x^, x) 



dL w dk p 
dk^ dx v 



+ 



dL w 
dx v 



dx v = -Vk{A) + 



dL w 
dx 







where Vk(A) = Va^ = A^V^k. Now let L w (k,x) be the dispersion relation 
on the 8-dimensional phase space (k,x), with both k and x as independent 
coordinates. The ray equations are 



dk 
da 



Vk{A) 



dL w 
dx 



and 



dx 
da 



Vx{A) — A — — 



dL w 
dk 



(43) 



where a is a parameter along the ray. Equation (J43|) is the 4D canonical 
Hamiltonian form of the ray equations for the waves. Since A = Ng we can 
divide the ray equations by iV = dL w /du =>■ 



Vk{g) 



1 dL w 
N dx 



du 
dx 



and Vx(g) = g 



1 dL w 
'N dk 



A du 
N = U+ d^ 



Its always useful to write the equations in coordinate free form, since then 
they can be evaluated in any CS. 
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11 The Energy Momentum Tensor (EMT) 



The total canonical EMT for a system consisting of a slowing varying back- 
ground flow interacting with waves is 

dL 

T = Lg — DL — k— where L = L (z\ dz 1 ; x) + L w (z\ dz\ a, k; x) 

OK 

The EMT T divides naturally into 3 parts: 

T = T + I + W = background + interaction + waves, where (44) 

8L 

T = L g-DL ; I = —DL W ; W = -k^ = kA 

T is the EMT of the background alone, since D = dz k £g> d/ddz k and L 
depend only on the background fields. Since I is the derivative of L w with 
respect to the background fields, it is the EMT of the wave-background 
interaction. W involves only the wave field k and L w , so it is the EMT of 
the waves alone. 

We need to know the action of D on all the quantities in L: 



Dn = 


nP 


DJ = 


Jg - g J 


Du = 


-Pu 


Duo = 


KU 


Dk = 


UlPu — UKU 


Dk = 


D(uju + k) = 


D(k ■ a) = 


uau 


DT = 


(7 - 1)TP 


DB = 


BP- PB- uBu 


DB 2 = 


2(B 2 P - BB) 


D(B-k) = 


(B ■ k)P - kB + loBu 


D(B-a) = 


(B-a)P-aB 


Df = 


fc 2 P 


£>7 = 


-CiP 


Dc 2 = 


2 7 'c 2 P 


c = 


T dc v 



(45) 



and i = (7 - 1 - c 2 - 0/2. 

11.1 The Proper Components of the EMT 

The decomposition of the parts of the EMT with respect to u yields the 
proper energy and momentum density and flux. The wave EMT W is 

W = k ® A = N(oou + k)®(u + v g ) = Euu + Euv g + Nku + Nkv 9 (46) 
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The total EMT due to the waves is T w = I + W = the wave + interaction 
EMT. Since L depends only on scalar fields, T, T , and T w are symmetric, 
but I and W are not. We have T = To + T w , where To is the EMT of the 
background given in Part 1, and 

T w = I + W = E(uu + uv g + v g u) + NnVg + PIP (47) 

11.1.1 The Components of the EMT for Waves 

The energy and momentum density and flux for the waves are the components 
of W and T w . They can be read off from equations (f4lj|) and (|47j). The wave 
energy density is 

W(u, u) = T w (u, u)= E = Nu 

since the interaction energy density I(u,u) = 0. The interaction energy flux 
density is P(I(u)) = also. The wave momentum density is — P{W\u)) = 
Nk, as one would expect. However, this term in W is canceled by one in I 
to give the total wave + interaction momentum density 

-P(Tl(u)) = -P(P(u) + W\u)) = Ev g = -P(T w (u)) 

as one would expect. The wave pressure tensor is 

PT W P = PIP + PWP = PIP + Nnv g 

which is always symmetric. 

11.2 Sound Waves 

The EMT of the background is T = puu +pP. For the waves, v p = v g = ck, 
and 

T w = I + W = E{uu + cuk + cku + kk + j'P) 
The components of T w in a CMCS are 

T w (u, u) = E = Nuj -P{T w {u)) = -P(Tl(u)) = Eck 

PT W P = E{j'P + kk) 

The components of the total EMT are 

T = T + T w = (p + E)uu + Ec[uk + ku) + Ekk + (p + i)P (48) 
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11.3 MHD Waves 

The EMT of the background for MHD was given in equation (|19|). 



11.3.1 Alfven Waves 

The phase and group velocities are v p = va cos 9k and v g = v A B (see eq. 
(J33J) and §10.2.1. Using gSJ) to calculate the EMT yields 

T w = E(uu + v A uB + v A Bu) + wV[n/(l - c 2 )P + 2BB]/A 
T = T + T W 

11.3.2 Magnetosonic Waves 

The phase velocity, group velocity and energy density of the fast and slow 
magnetosonic modes are given in §9.2.2 and §10.2.2. The EMT for the waves 
is 

T w = E(uu + uv g + v g u) + PT W P 

The expression for the wave pressure tensor PT W P is quite complex and I 
will not reproduce it here. 

11.4 The Equations for The Total Wave-Background 
System 

The following equations completely determine the waves and the background 
flow: 

V-T = V- J = V- A = 

and the ray equations (|41j) and (|43|) . One must include the EMT of the 
waves T w when it is comparable to T , i.e., when the wave energy density E 
is comparable to the thermal or magnetic field energy density p or B 2 , e.g. 
in the solar wind. Waves transport energy and momentum in ways in which 
other processes cannot. A similar treatment can be given for plasma waves 
(Jacques 1988) and wave in other continua. 
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